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Abstract 

We present a rigorous mathematical treatment of the zero-field orbital magnetic suscep- 
tibility of a non-interacting Bloch electron gas, at fixed temperature and density, for both 
metals and semiconductors/insulators. In particular, we obtain the Landau-Peierls formula 
in the low temperature and density limit as conjectured by T. Kjeldaas and W. Kohn in 1957. 
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1 Introduction and the main results 



Understanding the zero-field magnetic susceptibility of a Bloch electron gas is one of the oldest 
problems in quantum statistical mechanics. 

The story began in 1930 with a paper by L. Landau |30| . in which he computed the diamagnetic 
susceptibility of a free degenerate gas. (Note that the rigorous proof of Landau's formula for free 
electrons was given by Angelescu et al. [T] and came as late as 1975.) For Bloch electrons (which 
are subjected to a periodic background electric potential), the problem is much harder and -to our 
best knowledge- it has not been solved yet in its full generality 

The first important contribution to the periodic problem came in 1933, when R. Peierls [34) 
introduced his celebrated Peierls substitution and constructed an effective band Hamiltonian which 
permitted to reduce the problem to free electrons. Needless to say that working with only one 
energy band instead of the full magnetic Schrbdinger operator is an important simplification. 
Under the tight-binding approximation he claimed that the dominant contribution to the zero-field 
orbital susceptibility of a Bloch electron gas in metals (at zero temperature) is purely diamagnetic 
and is given by the so-called Landau-Peierls formula which consists of replacing in the Landau 
formula the mass by the effective mass of the electron. He showed as well the existence of another 
contribution which has no simple interpretation and whose magnitude and sign are uncertain. 

In 1953, E.N. Adams [5] claimed that the Landau-Peierls susceptibility is not always the 
dominant contribution to the zero-field orbital susceptibility. By considering the case of 'simple 
metals' (for which the tight-binding approximation is not appropriate), he showed that there exists 
others contributions (certain have even positive sign!) coming from the bands not containing the 
Fermi energy. Besides in special cases these contributions are of the same order of magnitude than 
the Landau-Peierls formula. However any general formulation of these contributions were stated. 

In 1957, T. Kjeldaas and W. Kohn [26] were probably the first ones who suggested that for 
'simple metals' the Landau-Peierls approximation is only valid in the limit of weak density of 
electrons, moreover, the Landau-Peierls formula (see below (|1.17[) and (|1.18p ) has to be corrected 
with some higher order terms in the particle density, and these terms must come from the bands 
not containing the Fermi energy. 

These three papers generated a lot of activity, where the goal was to write down an exact 
expression for the zero-field magnetic susceptibility of a Bloch electron gas in metals at zero 
temperature. In what follows we comment on some of the most important works. 

The first attempt to address the full quantum mechanical problem -even though the carriers 
were boltzons and not fermions- was made by J.E. Hebborn and E.H. Sondhcimer [21] [22] m I960. 
Unlike the previous authors, they developed a magnetic perturbation theory for the trace per unit 
volume defining the pressure. The biggest problem of their formalism is that they assumed that 
all Bloch energy bands are not overlapping (this is generically false; for a proof of the Bethe- 
Sommerfeld conjecture in dimension 3 see e.g. [23]), and that the Bloch basis is smooth in the 
quasi-momentum variables. This assumption can fail at the points where the energy bands cross 
each other. Not to mention that no convergence issues were addressed in any way. 

In 1962, L.M. Roth [36] developed a sort of magnetic pseudodiffcrential calculus starting from 
the ideas of Peierls, Kjeldaas and Kohn. She used this formalism in order to compute local traces 
and magnetic expansions. Similar results are obtain by E.I. Blunt 0. Their formal computations 
can most probably be made rigorous in the case of simple bands. 

In 1964, Hebborn et al. [20] simplified the formalism developed in [22] and gave for the first 
time a formula for the zero-field susceptibility of a boltzon gas. Even though the proofs lack any 
formal rigor, we believe that their derivation could be made rigorous for systems where the Bloch 
bands do not overlap. But this is generically not the case. 

The same year, G.H. Wannier and U.N. Upadhyaya [40] go back to the method advocated 
by Peierls, and replace the true magnetic Schrodinger operator with a (possibly infinite) number 
of bands modified with the Peierls phase factor. They claim that their result is equivalent with 
that one of Hebborn and Sondheimer [22], but no details are given. Anyhow, the result uses in 
an essential way the non-overlapping of Bloch bands. At the same time, L. Glasser [TH] gave an 
expression of the bulk zero-field susceptibility in terms of effective mass by the usual nearly free 
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electron approximation. 

In 1969, P.K. Misra and L.M. Roth [32] combined the method of [35] with the ideas of Wannier 
in order to include the core electrons in the computation. 

In 1972, P.K. Misra and L. Kleinman |31j had the very nice idea of using sum-rules in order 
to replace derivatives with respect to the quasi-momcntum variables, with matrix elements of the 
"true" momentum operator. They manage in this way to rewrite the formulas previously derived 
by Misra and Roth (which only made sense for non-overlapping bands) in a form which might also 
hold for overlapping bands. 

As we have already mentioned, the first serious mathematical approach on the zero-field sus- 
ceptibility appeared as late as 1975, due to Angelescu et al. PQ. Then in 1990, B. Hclffcr and J. 
Sjostrand [53] developed for the first time a rigorous theory based on the Pcicrls substitution and 
considered the connection with the dc Haas- Van Alphen effect. These and many more results were 
reviewed by G. Nenciu in 1991 [33]. A related problem in which the electron gas is confined by a 
trapping potential was considered by M. Combescure and D. Robert in 2001 [TT]. They obtained 
the Landau formula in the limit h — > 0. 

Finally we mention that the magnetic response can be described using the semiclassical theory 
of the orbital magnetism and the Berry-phase formula, see |35] for further details. The link between 
this approach and our work has yet to be clarified. 

Our current paper is based on what we call magnetic perturbation theory, as developed by the 
authors and their collaborators in a series of papers starting with 2000 (see [T3J [JJ] HS1 HI |3J 
|4j [5] [6] [7] and references therein) . The results we obtain in Theorem 11.21 give a complete answer 
to the problem of zero-field susceptibility. Let us now discuss the setting and properly formulate 
the mathematical problem. 

1.1 The setting 

Consider a confined quantum gas of charged particles obeying the Fermi-Dirac statistics. The 
spin is not considered since we are only interested in orbital magnetism. Assume that the gas is 
subjected to a constant magnetic field and an external periodic electric potential. The interactions 
between particles are neglected and the gas is at thermal equilibrium. 

The gas is trapped in a large cubic box, which is given by A^ = (— ^,^),L>1. 

Let us introduce our one-body Hamiltonian. We consider a uniform magnetic field B = (0, 0, B) 
with B > 0, parallel to the third direction of the canonical basis of R 3 . Let a(x) be the symmetric 
(transverse) gauge a(x) := 1(— x 2 , %i, 0) which generates the magnetic field (0,0,1). 

We consider that the background electric potential V is smooth, i.e. V <G C°°(IR 3 ) is a real- 
valued function and periodic with respect to a (Bravais) lattice T with unit cell £7. Without loss 
of generality, we assume that T is the cubic lattice Z 3 , thus 12 is the unit cube centered at the 
origin of coordinates. 

When the box is finite i.e. 1 < L < oo, the dynamics of each particle is determined by a 
Hamiltonian defined in L 2 (A^) with Dirichlct boundary conditions on dA^: 

H l {cj) = i( - iV x - c^a(x)) 2 + V L (x) (1.1) 

where Vl stands for the restriction of V to the box A^. Here lj := |£? G R denotes the cyclotron 
frequency. The operator Hl(uj) is self-adjoint on the domain D(Hl(u)) = Hq(Al) HH 2 (Al). It is 
well known (see [37]) that Hl{uj) is bounded from below and has compact resolvent. This implies 
that its spectrum is purely discrete with an accumulation point at infinity. We denote the set of 
eigenvalues (counting multiplicities and in increasing order) by {ej(w)}j>i. 

When L = oo we denote by Hoo(w) the unique self-adjoint extension of the operator 

i(^V x -wa(x)) 2 + y(x) (1.2) 

initially defined on C^(R 3 ). Then H^lo) is bounded from below and only has essential spectrum 
(see e.g. [8]). 
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Now let us define some quantum statistical quantities related to the quantum gas introduced 
above. For the moment we use the grand canonical formalism. The finite volume pressure and 
density of our quantum gas at inverse temperature /3 := (fc^T) -1 > (ks stands for the Boltzmann 
constant), at fugacity z := c"^ > ((i £ 1 stands for the chemical potential) and at cyclotron 
frequency uj G ffi are given by (see e.g. [H]): 

1 °° 

ft 03, z,u>) := In (1 + ze-P**M)} =^^(1 + ze~^) (1.3) 

Pi (/3,z, W ) :=Pz—{p,z, U ) = \——^- y (1.4) 

i=i 

As the semi-group &~^ Hl ^ is trace class, the series in (jl.3l) and (|1.4[) are absolutely convergent. 
Since the function K 9 oj n- P L (f3,z,u!) is smooth (see [7]), we can define the finite volume 
orbital susceptibility as the second derivative of the pressure with respect to the intensity B of 
the magnetic field at B = (see e.g. [T]): 

Xj? c (J3,z):= (l) ^(P> z >°)- (1-5) 

When fills the whole space, we proved in [33] that the thermodynamic limits of the three 
grand canonical quantities defined above exist. By denoting P oc (/3, z,o;) := limL_ s . 00 Pl{/3, z, uj), 
we proved moreover the following pointwise convergence: 

d Poo d Pl 

poo(/3, z,u) := (3z— —(P,z,w) = lim (3z— — (/?, z,w) (1.6) 

OZ L— >oo OZ 

and the limit commutes with the first derivative (resp. the second derivative) of the grand canonical 
pressure with respect to the fugacity z (resp. to the external magnetic field B). 

Now assume that our fixed external parameter is the density of particles po > 0. We prefer to 
see poo as a function of the chemical potential p instead of the fugacity z\ the density is a strictly 
increasing function with respect to both p and z. Denote by Poo{P,po) G K the unique solution 
of the equation: 

Po=Poo(/3,e^~^ o) ,o). (1.8) 
The bulk orbital susceptibility at (3 > and fixed density pa > defined from (11.71) is defined as: 

X{f3,p Q ) := Xg c (p,e^°°^) . (1.9) 

In fact one can also show that X(f3,p ) = — (f) 2 ^? 3 (P,po,0) where foo(P, Po, w) is the thermo- 
dynamic limit of the reduced free energy defined as the Legendre transform of the thermodynamic 
limit of the pressure (see e.g. [38]). Note that for a perfect quantum gas and in the limit of low 
temperatures, (|1.9[) leads to the so-called Landau diamagnetic susceptibility, see e.g. [1]. 

In order to formulate our main result, we need to introduce some more notation. In the case 
in which w = 0, the Floquet theory for periodic operators (see e.g. [TU], [55] and section 3) allows 
one to use the band structure of the spectrum of Hoo(0). Denote by fl* = 27rf2 the Brillouin zone 
of the dual lattice T* = 2nZ 3 . 

If j > 1, the jth Bloch band function is defined by £j := [minkgn* EjQc), maxkgu* Ej(\s)] 
where {Ej(k)}j>i is the set of eigenvalues (counting multiplicities and in increasing order) of the 
fiber Hamiltonian h(k) := |(-iV + k) 2 + V living in L 2 (T 3 ) with T 3 := M 3 /Z 3 the 3-dimensional 
torus. With this definition, the Bloch energies Ej(-) are continuous on the whole of 12*, but they 
are differentiable only outside a zero Lebesgue measure subset of 12* corresponding to cross-points. 
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In the following we make the assumption that the Ej 's are simple eigenvalues for k in a subset of 
n* with full measure. Note that this assumption is not essential for our approach but it simplifies 
the presentation, sec Remark 3 below the Theorem 11.21 

The spectrum of -ffoo(O) is absolutely continuous and given (as a set of points) by <7(-ffoo(0)) = 
U^Li £j- Note that the sets £j can overlap each other in many ways, and some of them can even 
coincide even though they are images of increasingly ordered functions. The energy bands are 
disjoint unions of £j's. Moreover, if max£j < min^j+i for some j > 1 then we have a spectral 
gap. Since the Bethe-Sommcrfcld conjecture holds true under our conditions [23], the number of 
spectral gaps is finite, if not zero. 

It remains to introduce the integrated density of states of the operator if^O). Recall its 
definition. For any let Nl(E) be the number of eigenvalues of Hl(0) not greater than E. 

The integrated density of states of -ffoo(O) is defined by the limit (see [TT]): 

r ™ r Mg> i- Tr {x(-oo,£](gL(0))} 
noo(-E) := km rr , = 1™ — nrn ( L1 °) 

L-s-oo |Ai| L-s-oo | Ai| 

and n oc (-) is a positive continuous and non-decreasing function (see e.g. |10j). In this case one 
can express (E) with the help of the Bloch energies in the following way: 

noo{E) = ^^j^ XM (^.(k)) dk (1.11) 

where X[_b ,b](') ^ s the characteristic function of the interval [Eq, E]. Thus is clearly continuous 
in E due to the continuity of the Bloch bands. Moreover, this function is piecewise constant when 
E belongs to a spectral gap. 

1.2 The statements of our main results 

The first theorem is not directly related to the magnetic problem, and it deals with the rigorous 
definition of the Fermi energy for Bloch electrons. Even though these results are part of the 'physics 
folklore', we have not found a serious mathematical treatment in the literature. 

Theorem 1.1. Let pq > be fixed. If fioo{PiPo) * s the unique real solution of the equation 
Poo (/3,e' 3 '\0) = po (see (|1.8[) ). then the limit: 

£f(po) ■= lim /!oo(/3,Po) (1-12) 

13— too 

exists and defines an increasing function of po called the Fermi energy. There can only occur two 
cases: 

SC (semiconductor /insulator /semimetal): Suppose that there exists some N G N* such that 
Po = n aD (E) for all E £ [max£/v, minfjv+i]- Then: 

maxfjv + mm£ N+1 
£f(Po) = 2 • ( 3 ) 

M (metal): Suppose that there exists a unique solution Em of the equation n ao (EM) = Po 
which belongs to (minfjy, max£jv) for some (possibly not unique) N. Then : 

£f{po) = E M . (1-14) 



Remark 1 . In other words, a semiconductor/semimetal either has its Fermi energy in the middle 
of a non-trivial gap (this occurs if max£/v < min£/v + i), or where the two consecutive Bloch bands 
touch each other closing the gap (this occurs if max£/v = rnin£/v + i). As for a metal, its Fermi 
energy lies in the interior of a Bloch band. 
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Remark 2. According to the above result, £p is discontinuous at all values of po for which the 
equation n oa (E) = po does not have a unique solution. Each open gap gives such a discontinuity. 

Now here is our main result concerning the orbital susceptibility of a Bloch electrons gas at 
fixed density and zero temperature: 

Theorem 1.2. Denote by E := inf (7(^(0)) . 

(i). Assume that the Fermi energy is in the middle of a non-trivial gap (see (|1.13[) )■ Then 
there exist 2N functions Cj(-),T)j(-), with 1 < j < N, defined on f2* outside a set of Lebesgue 
measure zero, such that the integrand in (|1.15p can be extended by continuity to the whole of Q* 
and: 



p— >oo 



2 -/V 

^) 5(2^ ^dk^lc^ + ^OcJ-MPo)}^-^)}. (1.15) 



(ii). Suppose that there exists a unique N > 1 suc/i i/iai £f{po) € (min£/v, max£/v). Assume 
that the Fermi surface Sf '■= {k G 0* : -E/v(k) = £f(/9o)} is smooth and non-degenerate. Then 
there exist 2N + 1 functions J r /v(-), Cj (•), t)j (■) uraf/i 1 < j 1 < iV, defined on Q,* outside a set 
of Lebesgue measure zero, in such a way that they are all continuous on Sf while the second 
integrand in (|1.16[) can be extended by continuity to the whole of f2* ; 



Xm{po) 



lim X(J3 tPo ) = - 



1 1 



c) 12 (2tt) 3 



(l.K 



- 6 



dg(k) 
Sf \VE N (k) 

N 

dkE 



a 2 £;Ar(k) d 2 E N (k) fd 2 E N (k)\ 



dk 2 



dk\dk2 I 



- 3-FAr(k) 



X[ Eo ,£ F ( PO )] (Ej(k))cj(k) + {Ejik) - £ F (pn)}x[E ,£ F (po)]( E j( k ))^j(. k ) 



Here X[e .£f(po)](' ) denotes the characteristic function of the interval Eq < t < £p(po)- 

(iii). Let kp := (6ir 2 po)i be the Fermi wave vector. Then in the limit of small densities, (| 1 . 1 6j) 
gives the Landau- Peierls formula: 



Xm(po) 



24ir 2 c 2 



o(k F )\ 



(1.17) 



[?] i<i<3 are ^ e e ^9 enva ^ ues of the positive definite Hessian matrix {9 I 2 :) -£ , i(0)}i<ij<3. 



Remark 1 . The functions Cj(-) and with 1 < j < N which appear in (|1.15[) are the same as 
the ones in (|1.16[) . All of them (as well as .Fjv(-)) can be explicitely written down in terms of Bloch 
energy functions and their associated eigenfunctions. One can notice in (|1.16[) the appearance of 
an explicit term associated with the Nth Bloch energy function; it is only this term which will 
generate the linear kp behavior in the Landau-Pcicrls formula. 

Remark 2. The functions Cj(-) and Oj(-) might have local singularities at a set of Lebesgue 
measure zero where the Bloch bands might touch each other. But their combinations entering the 
integrands above arc always bounded because the individual singularities get canceled by the sum. 
Remark 3. The results in (i) and (ii) hold true even if some Bloch bands are degenerate on a 
subset of full Lebesgue measure of fl*. But in this case the functions Cj(-), Dj(-) and Tn(- ) cannot 
be expressed in the same way as mentioned in Remark 1. Their expressions are more complicated 
and require the use of the orthogonal projection corresponding to Ej{-), see the proof of Lemma 
13.71 for further details. 
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Remark 4 ■ When m* = m\ — m\ = m* holds in (iii) , (|1.17|) is nothing but the usual Landau- 
Pcicrls susceptibility formula: 



X M (po) ~ -7TT- o — r^ k F when k F -> 0. (1-18) 

Note that our expression is twice smaller than the one in |34j since we do not take into account 
the degeneracy related to the spin of the Bloch electrons. 

Remark 5. The assumption V e C°°(T 3 ) can be relaxed to V e C r (T 3 ) with r > 23. The 
smoothness of V plays an important role in the absolute convergence of the series defining X(J3, po) 
in Theorem 13. 11 before the zero-temperature limit; see |16) for a detailed discussion on sum rules 
and local traces for periodic operators. 

Remark 6. The role of magnetic perturbation theory (see Section 3) is crucial when one wants 
to write down a formula for X(/3,po) which contains no derivatives with respect to the quasi- 
momcntum k. Remember that the Bloch energies ordered in increasing order and their corre- 
sponding eigenfunctions arc not necessarily differentiable at crossing points. 

Remark 7. We do not treat the semi-metal case, in which the Fermi energy equals £f(pq) = 
max£/v = min£/v+i for some N > 1 (see (|1.13[l ). This remains as a challenging open problem. 



1.3 The content of the paper 

Let us briefly discuss the content of the rest of this paper: 

• In Section 2 we thoroughly analyze the behavior of the chemical potential poc when the 
temperature goes to zero defining the Fermi energy. These results are important for our 
main theorem. 

• In Section 3 we give the most important technical result. Applying the magnetic perturbation 
theory we arrive at a general formula for X(/3, po) which contains no derivatives with respect 
to k. The strategy is somehow similar to the one used in [14] for the Faraday effect. 

• In Section 4 we perform the zero temperature limit and separately analyze the situations in 
which the Fermi energy is either in an open spectral gap or inside the spectrum. It contains 
the proofs of Theorem II .21 (i) and (ii). 

• In Section 5 we obtain the Landau- Peier Is formula by taking the low density limit. It contains 
the proof of Theorem ll.2l (iii). 



2 The Fermi energy 

This section, which can be read independently of the rest of the paper, is only concerned with 
the location of the Fermi energy when the intensity of the magnetic field is zero (i.e. ui = 0). In 
particular, we prove Theorem 11.11 Although we assumed in the introduction that V € C°°(T 3 ), 
all results of this section can be extended to the case V £ L°°(T 3 ). 



2.1 Some preparatory results 

Let £ >-¥ f(/3, p; £) := In (l + c^ - ^) be a holomorphic function on the domain {£ 6 C : 9£ G 
(—77/(8, n/P)}. Let r the positively oriented simple contour included in the above domain defined 
by: 

(2.1) 



r := {k£ e [5,oo), 3£ = ±^} u = 5,$*te 



7T 7T 

2/3' 2/3 



}• 



where S is any real number smaller than Eq := inf o^iJo^O)) < inf c(.ffoo(u;)). In the following we 
use S := Eq — 1. 
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The thermodynamic limit of the grand-canonical density at (3 > 0, p € K and to > is given 
by (see e.g. [8]): 

Poo^ J e^w) = |^^^ (R , ) |xn^f J , D 08 J /i;O(^ooH-O" 1 } ( 2 - 2 ) 

where )\f_d(/3, P\ : = — P~ X 9ff(f3, P\ £) ~ (c' 3 ^"^' + is the Fermi-Dirac distribution function 
and xo denotes the characteristic function of 17. We prove in (35] (even for singular potentials) 
that Poo(/3, •, w) can be analytically extended to the domain C \ (— oo, —e^ E °^]. 

Now assume that the intensity of the magnetic field is zero (uj = 0). The following proposition 
(stated without proof since the result is well known), allows us to rewrite (|2.2[) only using the 
Bloch energy functions k i— ► Ej(\t) of i?oo(0): 

Proposition 2.1. Let (3 > and /igl. Denote by fi* the first Brillouin zone of the dual lattice 
2ttZ 3 . Then: 

1 °° f 

Poo (/3,e^,0) = -— 3^ / dkfroGS, ,*;£,•(!£)). (2.3) 



Note that another useful way to express the grand-canonical density at zero magnetic field 
consists in bringing into play the integrated density of states (IDS) of the operator -ffoo(O) (see 
([TTTU|> for its definition): 

Poo{P,e^,Q) = - j°° dA^G9,wAKo(A). (2.4) 

When the density of particles po > becomes the fixed parameter, the relation between the 
fugacity and density can be inverted. This is possible since for all f3 > 0, the map Poc(/3, - ,0) is 
strictly increasing on (0,oo) and defines a C°°-diffeomorphism of this interval onto itself. Then 
there exists an unique Zoo{P,po) € (0, oo) and therefore an unique fj,oo(f3,po) € R satisfying: 

Poo (/3,e^M,0)=p . (2.5) 

We now are interested in the zero temperature limit. The following proposition (again stated 
without proof) is a well known, straightforward consequence of the continuity of noo(-): 



Proposition 2.2. Let p > Eq := inf o^iJo^O)) be fixed. We have the identity 
lim Poo ((3, e ^,0) = - 

p— >oo 12 



1 f 

lim Poo (/?,e^,0) = 77 —Y J <fcX[E ,»](E j (k))=n 0o (p), (2.6) 



j 

where X[E ,fi](') denotes the characteristic function of the interval [Eo,p] 



We end this paragraph with another preparatory result concerning the behavior of near 
the edges of a spectral gap. This result is contained in the following lemma: 

Lemma 2.3. Let P q > be fixed. Assume that there exists N > 1 such that n 00 (£') = po for all E 
satisfying max£/v < E < min£/v + i. We set aj/ := max£/v and bjy := rnin£/v+i- Assume that the 
gap is open, i.e. < &j\r- Then for S > sufficiently small, there exists a constant C = Cs > 
such that: 

n oo{ a N) — n ooW > C(ajv — A) 3 whenever A e [ajy — S, ajv] (2.7) 

and 

noo(A) — ?1oo(&a0 > C(A — &at) 3 whenever A € [bjy, £>jv + <$]. (2.8) 
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Proof. We only prove (12.71) . since the proof of the other inequality (|2.8|) is similar. Since ajy = 
maxkgn* -E/v(k), the maximum is attained in a (possibly not unique) point ko, i.e. = -Ejv(ko). 
This means that a at is a discrete eigenvalue of finite multiplicity 1 < M < N of the fiber operator 
h(ko) = h(— iV + ko) 2 + V. In particular, a^r is isolated from the rest of the spectrum since we 
assumed that ajy < &jv < -Ejv+i(ko). Now when k slightly varies around ko, the eigenvalue ajv 
will split into at most M different eigenvalues, the largest of which being i?Ar(k). Thus from the 
second equality in (|2.6p we obtain: 

noo(a N ) - rioo(A) > — !-^Vol{k € fT : A < i?jv (k) < 0^}. 

We now choose 8 small enough such that 

a{h{k )) n [a N -S,a N + S] = {E N (ko)}. 

We use analytic perturbation theory in order to control the location of the spectrum of h(k) when 
|k — ko | is small (we assume without loss of generality that ko lies in the interior of Q*). By 
writing 

h(k) = ftfko) + (k - k ) ■ (-iV + k ) + (k - k ) 2 /2 =: h(k ) + W(k), 
we see that we can find a constant C > such that 

\\w(k)(h(k ) - i)- 1 !! < c|k - k |, |k - k | < i. 

Take a circle 7 with center at a at and radius r = (a^ — A)/2 < (5/2. For any z £ 7, by virtue of 
the first resolvent equation: 

(h(k ) - z)" 1 = (/i(ko) - i)" 1 + 0* - i)(/i(k ) - *)- 1 (^(k ) - z)- 1 

and by using the estimate ||(/i(ko) — ) 1 1 1 = 2/{a^ — A), we can find another constant Cg > 
such that: 

sup \\w{k)(h(k Q ) - z)" 1 !! < c 5 j k ~ k °' |k - ko| < 1. 

It turns out that if |k — ko|/(aAr — A) is smaller than some e > 0, then 

sup \\W(k)(h(k ) - z) _1 || < eC s whenever |k- k | < {a N - A)e. 

Standard analytic perturbation theory insures that if e is chosen small enough, h(k) will have 
exactly M eigenvalues inside 7. Thus for all k satisfying |k — k | < e(a.N — A), we have a(h(k)) fl 
[a^ — S, ajv] C [(aw + A)/2,aAr] C [A,aAr]. In particular, A < E^v(k) < for all such k's. But 
the ball in Q* where |k — k | < (a» — A)e has a volume which goes like (a^ — A) 3 , and the proof 
is finished. 



2.2 Proof of Theorem [TTT1 

In this paragraph we prove the existence of the Fermi energy. We separately investigate the 
semiconducting case and the metallic case. 

2.2.1 The semiconducting case (SC) 

We here consider the same situation as in Lemma 12.31 in which there exists N > 1 such 
that n ao (E) = po for all E satisfying max£/v < E < minf/v+i- We set ajy := max£/v and 
6at := min^Ar+i. Let p(f3) := Poo{P, Po) be the unique solution of the equation poo(fl, e^*, 0) = po- 
We start with the following lemma: 
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Lemma 2.4. 



o-n < Pi := liminf p(ft) < limsup/i(/3) =: ji2 < &jv- (2-9) 

P^oo /3->oo 



Proof . We will only prove the inequality < fix, since the proof of the other one (fi 2 < &/v) is 
similar. Assume the contrary: /ii < ajv- Define e := — /^i > 0. Then there exists a sequence 
{(3n}n>i with ft n — > oo and an integer M e > 1 large enough such that: 

lim p(ft n ) = /ii and p(ft n ) < — e/2 < ajy, Vn > M e . 

n— >oo 

Since p oc (/3, e' 3 ^, 0) is an increasing function of p, we have: 

By letting n — > oo in the above inequality, (|2.6j) implies: 

Po < n-oo(ajv - e/2) < noo(ajv) = p 
where in the second inequality we used (|2.7[) . We have arrived at a contradiction. ■ 



Now if ajv = frjv, the proof of (|1.13p is over. Thus we can assume that ajy < 6jv, i.e. the gap 
is open. We have the following lemma: 

Lemma 2.5. Define cjy = (on + &a0/2. For any < e < (6jy — ajv)/2, £/iere exists ft e > Zarge 
enough such that p(ft) G [cjy — e, c/v + e] whenever ft > ft e . 

Proof . We know that p(ft) exists and is unique, thus if we can construct such a solution in the 
given interval, it means that this is the one. We use (|2.4p in which we introduce p(ft) and arrive 
at the following identities: 

dX (ft, p(ft); X)n (Xl (X) = - / dX (ft, /i(/3); A)n oc (A) 
-oo 9X J_ x <9A 

- «oo(&w)fF_D(/3,A i (/3); +«oo(aAr)fFD(/3,^(/3);aAr) 

where in the last term we used the fact that rioo(') is constant on the interval [aiy,b]\r], and this 
constant is nothing but po- We can rewrite the above equation as: 

dX^^(ft,p(ft);X){n 00 (Xyn 00 (a N )} ^ f dX^^(ft, p(ft); X){ noo (b N )- noo (X)} (2.10) 
<9A J bN dX 

where we used the fact that fFij(/3, p(ft)\ — oo) = 1 and fFD(ft, p(ft)', A) < Ce~ x/3 for large A. 
In the left hand side of (|2.10[) we now introduce the explicit formula: 

d x f FD (ft,p(ft);X) = -P {Q ^ m + l)2 = -^ (aN -» m {e ^ m + ir 
while in the right hand side of (|2.10j) we use another expression: 

-/3(A- M (/3)) e -/3(A-b„) 

9\f F D(ft, m(/3); A) = -ft——-—— = -fr-fito-Hn) 



(1 + C -/3(A- M (/J)))2 ^ (l + e -«M«)) 2 ' 

Then (|2.10p can be rewritten as: 

°™ e /3(A-a N ) 

eMwMa ™ +bN)} y dA ( i +c - g (A-, W ))2 ^( A ) - »«m>. ( 2 - n ) 
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or by taking the logarithm: 



If ( e /3(A-ajv) 
r°° e -/3(A-6 N ) 

Let us define the smooth function / : [cat — e, cat + e] 1— > R given by: 



ln( / rfA ,,, , „_^ A _„.rm^ 2 { n °°( A )~ n °°( 6j v)> 



/(*) := c N + — I In ^ y ^ rf A (c/j(A _ :c) + i)2 {n 00 (a j v) - n x (X)} 
( r°° P -/3(A-b„) \ ~| 

- ln {J bN dX -^T^^ Mx) ~ noo(bN)} )\- (2 ' 12) 

We will prove that if ft is large enough, then / invariates the interval [cat — e, cat + e], which is 
already enough for the existence of a fixed point. This would also show that p(ft) must be in that 
interval. But in fact one can prove more: / is a contraction for large enough ft. 

The idea is to find some good upper and lower bounds when ft is large for the integrals under 
the logarithms. We start by finding a lower bound in ft for the first integral. Let S > sufficiently 
small. Using (|2.7[) in the left hand side of (|2.11[) we get: 

ajv p/3(A— ajv) f< f»» 

^{n 00 (a JV )-n 00 (A)}>- / e-^-^(a N -X) 3 (2.13) 

where we used that x > > A in order to get rid of the numerator. After a change of variables 
and using some basic estimates one arrives at another constant C > such that for ft sufficiently 
large: 

(•a N g/3(A-ajv) Q 

dX ( c g(A-x) + 1)2 i n «>M - n oc (A)} > . (2.14) 



00 



By restricting the interval of integration to [6 at, £>at + <5] and by using ()2.8[) . we obtain by the same 
method a similar lower bound for the second integral under the logarithm. Moreover, using the 
Weyl asymptotics which says that noo(A) ~ A^ for large A (see e.g. [27]), one can also get a 
power-like upper bound in ft for our two integrals. 

We deduce from these estimates that there exists a constant C £ > such that: 

sup \f(x)- CN \<^^., ft>l. 

x£[c N — e,c N +e] P 

Thus if ft is large enough, / invariates the interval. Being continuous, it must have a fixed point. 
Moreover, the derivative f'(x) decays exponentially with ft uniformly in x € [cat — e, cjv + e]. It 
implies that if ft is large enough, then ||/'||oo < 1, that is / is a contraction. ■ 



2.2.2 The metallic case (M) 

Consider the situation in which there exists a unique solution Em of the equation n 0o (-E'Jw) = 
Po, and this solution lies in the interior of a Bloch band. In other words, there exists (a possibly 
not unique) integer N > 1 such that min^ < Em < max^Ar. Wc will use in the following that 
the IDS Uoo(-) is a strictly increasing function on the interval [min£ at, maxfjv]. 

Let fi(ft) := Hoo(P,po) be the unique real solution of the equation Poc(ft, e^^\ 0) = po- Let 
us show that : 

E M < Uminf p,(ft) < limsup^(/3) < E M , (2.15) 
which would end the proof. We start with the first inequality. 
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Assume ad-absurdum that p\ := liminf^^oo p(/3) < Em- Then there exists e > and a 
sequence {/3„}„>i satisfying f3 n — ¥ oo such that: 

lim p(/3 n ) = p lt /Lt(/3„) < E M - e, Vn > 1. 

Since poc(/3, e^, 0) is increasing with we have: 

n^Eu) = Po = lim P ooWn,e MM ,0) < lim Poo (/3„, e A.< B "-0 0) = noo {E M - e), (2.16) 

where in the last equality we used (|2.6[) . But the inequality u^Em) < n oc (EM — e) is in 
contradiction with the fact that n QO (-) is a strictly increasing function near Em- Thus Em < 

Now assume ad-absurdum that p2 '■= lim sup^^ p{j3) > -Em- Then there exists e > and a 
sequence {/3 n } n >i satisfying (3 n — > oo such that: 

lim p(/3 n ) = p 2l E M + e < p(/3 n ), Vn > 1. 

n— yoo 

We again use that p oc (/3, e' 9 ' 1 , 0) is increasing with p and write: 

nooiEM + e) = lim Poo (/3„, e ^( s «+ £ ), 0) < lim Poo (fi n , e ^^\ 0) = p = ««, (2.17) 

where in the first equality we again used (|2.6p . But the inequality u^Em + e) < n oc (EM) is 
also in contradiction with the fact that n oc (-) is a strictly increasing function near Em- Therefore 
P2 < E M - ■ 

3 The zero-field susceptibility at fixed density and positive 
temperature 

In this section we prove a general formula for the zero-field grand-canonical susceptibility of a 
Bloch electrons gas at fixed density and positive temperature. 
Here is the main result of this section: 

Theorem 3.1. Let (3 > and po > be fixed. Let p^ — Poo{(3,po) € ^ the unique solution of 
the equation /^(/S, e^ ,u> = 0) = po- Then for each integer j\ > 1 there exists four families of 
functions ), with I G {0, 1,2,3}, defined on CI* outside a set of Lebesgue measure zero, such 

that the integrand below can be extended by continuity to the whole of CI* : 



X(P,po) 



\ 2 i i 
e \ 1 1 



33 E / dkj^Ui^^E^k^Ak), (3.1) 

I J1= l Jn * 1=0 



cj 2^(27 

with the convention (<9°f)(/3, p oc ;E J1 (k)) = f(/3, p oc ;Ej 1 (k)) := ln(l + e ^^~-- E «( k ))). 



This formula is a necessary step in the proof of Theorem II .21 (i) and (ii) (this is the aim of the 
following section) when we will take the limit of zero temperature. 

The special feature of this formula lies in the fact that each function Cj l! i(-) can be only 
expressed in terms of Bloch energy functions and their associated eigenfunctions. For each integer 
ji > 1, the functions Cj lt 2(-) and Cj 1 ,3(-) are identified respectively in p.28[) and (|3.27[) . As for the 
functions Cj lt i(-) with I G {0, 1}, they can also be written down but their explicit expression is not 
important for the proof of Theorem 11.21 Note as well that the above formula brings into play the 
Fermi-Dirac distribution and its partial derivatives up to the second order. This will turn out to 
be very important when we will take the limit (3 — > oo in the following section. 
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3.1 Starting the proof: a general formula from the magnetic perturba- 
tion theory 

We start by giving a useful formula for the thermodynamic limit of the grand-canonical sus- 
ceptibility. Let (3 > and z := e^ £ (0, oo) the fixed external parameters. Let L be the positively 
oriented contour defined in (|2.1[) , going round the half-line [E , oo), and included in the analiticity 
domain of the map £ i-> f(/3, z;£) = ln(l + ze~ 13 ^). Denote by i?oo(w,£) := {H^^) — £) _1 for all 
£ G /9(iJoo(w)) and to € R. Taking into account the periodic structure of our system, it is proved 
(see [7], Theorem 3.8) that the thermodynamic limit of the grand-canonical pressure of the Bloch 
electron gas at any intensity of the magnetic field B is given by: 

PooOM.w) == flL^-Tr i2(K3) (xn / dffC8,z;0«oo(w,o}. ( 3 - 2 ) 



/3|0| 2tt 



r 



where f2 is the unit cube centered at the origin of coordinates (xn denotes its characteristic 
function). Although the integral kernel i?oo(- , ■ ;^ 7 £) of the resolvent has a singularity on the 
diagonal, the integration with respect to £ in (|3.2[) provides us with a jointly continuous kernel on 
R 3 x R 3 . One can see this by performing an integration by parts in (|3.2[) and using the fact that the 
kernel of R^(u, £) is jointly continuous. Moreover, one can prove [H[5j[6] that the thermodynamic 
limit of the grand-canonical pressure is jointly smooth on (z,ui) € (— e^ E °,oo) x R. 

Let u G R and £ € p(H 00 (ui)). Introduce the bounded operators T 00 i(a;,£) and T 00i 2(u;,£) 
generated by the following integral kernels: 

Too,i( x >y; w >£) :=a ( x -y) ■ (*V x + wa(x)).R 00 (x,y;u;,£) (3.3) 

Too,2(x,y;u;,0 := ^a 2 (x - y)i? oc (x, y; £), x ^ y, (3.4) 

where a(-) stands for the usual symmetric gauge a(x) = A x = ^(—X2,xi,0). We introduce 
the following operators : 

Woo,iG8,^w) := / def(/3, W 0floo(w,0roo,i(w,0roo,i(w,0 (3.5) 
Woo,aG8,/*,w) ~ 7r~ [ d£f(ftw0 i *»(w>0r O o,3(w,f) (3.6) 

27T J r 

One can prove using the same techniques as in [15] that these operators are locally trace class and 
have a jointly continuous kernel on R 3 x R 3 . By a closely related method as the one in [4], [5], it 
is proved in |39j that we can invert the thermodynamic limit with the partial derivatives w.r.t. to 
of the grand-canonical pressure. Then the bulk orbital susceptibility reads as: 

^1^1 {Tr L 2 (K3) {xf2yVoo4(/3 I A i ,^)} - Tr L 2 (R3) {xoW o,2(/3,A 1 I w )}} 

We mention that the above formula is obtained using the so-called 'gauge invariant magnetic 
perturbation theory' applied to the resolvent integral kernel (see e.g. [T5] for further details) 
which allows to control the linear growth of the magnetic vector potential. 

The quantity which we are interested in is the orbital susceptibility at zero magnetic field and 
at fixed density of particles pq. Note that the pressure is an even function of u, thus its first order 
derivative at lo = is zero. This explains why the susceptibility is the relevant physical quantity 
for the weak magnetic field regime. 

The orbital susceptibility at zero magnetic field and fixed density po is given by (see also (|1.8|0 : 

fe\ 2 2 r 1 (3-7) 

= I -J ■g|^|Tr L 2 (R 3 ) {xnW O o4(/3:Moo,0)} - Tr L 2 (R 3 ) {xnW oo ,2(/?,A i oo,0)}|. 



e\ 2 2 



13 



The formula Q3.7[) constitutes the starting-point in obtaining (|3.1j) . The next step consists in 
rewriting the local traces appearing in (|3.7[) in a more convenient way: 

Proposition 3.2. Letp a := —id a with a £ {1, 2, 3} be the cartesian components of the momentum 
operator defined in L 2 (M. 3 ). Then we have: 



Tr L 2 (R 3 ) {xoW oo ,i(^,^oo,0)} = |2^Tr L 2 (R 3) jxn J d£f(/3, ^; £) 

[iZc (0, OpiiZoo (0, OPlRoc (0, f ) {p 2 i?oo (0, OPlRco (0,0- Pi i?oo (0, OP2R00 (o, } + 

+ R^iO, e)p 2 i?oo(0, OPi^co(0, 0{pii?co(0, 0P2i?oo(0, - p 2 i?co(0, Opi^co(0, £)}] } (3-8) 

and 

Tr^^lxnWoo^^j/iocO)} = -i^-Tr^s) jxn J d^f(f3, fi^;^) 

i?oo (0, £)i?oo (0, [p 2 i?oo (0, C)p 2 i?oo (0,6+ Pli?oo (0, OPlRoo (0,0- i?oo (0, 0] l ■ (3-9) 

Proof. We begin with the justification of (|3.9[) . By rewriting (|3.4p as: 

Too )2 (x, y; 0, £) = ^{e 3 A (x - y)} • {e 3 A (x - y)} J R co (x, y; 0, 
= g [(a* - y 2 ) 2 + (xi - yi) 2 ] Roo(x, y; 0, £), 

from p.6p it follows: 

Woo, 2 (x,x;/3 jA i,0) (3.10) 



i / d£f(J3,n;£) [ dzi 1 , 00 (x, Z ;0,0[(^ 2 -2: 2 ) 2 + (^i-a;i) 2 ]^oo(z,x;0,0, Vx e F 3 
° Jr Jr 3 



Let / £ {1,2}. Denote by X the multiplication operator with x. Then for all z 7^ x we can write: 

(^-x ( )-Roo(z,x;0,0 = [X.e ! ,.R oo (0 I 0](s,x) = {^(0,0 [#oo(0),X • e^M^x). 
We know that [iJoo(O), X • e;] = — ip;. Thus: 

(2 ; -x i )-R 00 (z,x;0,0 = -H ii °o(0,0^cx>(0,e)}(z,x). (3.11) 
Using standard commutation rules, we deduce from (|3.11[) that for I G {1, 2} and for all z/x: 
{zi - xifRoo (z, x; 0, £) = -{2Roo (0, f Jpi^oo (0, f )Pi^ (0, - #00 (0, £)i?oo (0, £) } (z, x). (3.12) 



It remains to put (|3.12| in (|3. 10|) . and we get (|3.9|) . 

Let us now prove now (|3.8|) . Since the divergence of a is zero, then for x ^ y we have: 



Too,! (x, y; 0, = \ V x • {e 3 A (x - y)}i? oc (x, y; 0, 

(x2 - 2/2) . {x\ - vi) 
7. e i H ^ e 2 



From (13.51) it follows that for all x G 



i?oo(x,y;0,0- 



W OOjl (x,x; ( S,/i,0) = - / d^f(/3,/i;0 / dz x / dz 2 R^x, z l5 0, 
4 Jr Jr 3 Jk 3 

{(«V Zl • ei)[-(zi ;2 - 2 2j2 )i? 00 (0,0(zi,z 2 )] + (iV ai • e 2 )[(z M - z 2 ,i)i?oo(0, 0( z i- z 2)] }' 

• {(«V Z2 ■ ei)[-(« 2 ,2 - a; 2 )^oo(0,0( z 2,x)] + (iV Z2 • e 2 )[(z 2ji - x 1 )R oo (0,^)(z 2 ,x)}}. 
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Then by using (|3.11[) . we get (|3.8[) from the following identity: 

Vxel 3 , W OOll (x,x;/3,/i,0) = i /d^fGS.wO / dz x / dz 2 i? oc (x, z i; 0, 

4 Jr it 3 Jr3 

{jpi(i?oo(0,C)P2i?oo(0,C))(zi,Z 2 ) - jp 2 (i?oo(0,OPl-Roo(0,0)(zi,Z2)} 

■ 

3.2 Using the Bloch decomposition 

We know that (see e.g. [10]) -ffoo(O) can be seen as a direct integral J®, dk/i(k) where the fiber 
Hamiltonians h(k) acting in L 2 (T 3 ) are given by : 

h(k) = i(-iV + k) 2 +V. (3.13) 

Recall that h(k) is essentially self-adjoint in C°°(T 3 ); the domain of its closure is the Sobolev 
space "H 2 (T 3 ). For each k£!l*, h(k) has purely discrete spectrum. We have already denoted by 
{Ej(h)}j>i the set of eigenvalues counting multiplicities and in increasing order. The correspond- 
ing eigenfunctions {uj(- ; k)}j>i form a complete orthonormal system in L 2 (T 3 ) and satisfy: 

/i(k)u,-(-;k) = Ej(lL)uj(-;lL). 

The eigenfunctions Uj 's are defined up to an arbitrary phase depending on k. These phases cannot 
be always chosen to be continuous at crossing points, and even less diffcrentiable. For the following 
let us introduce another notation. Let a £ {1, 2, 3}, and let i, j > 1 be any natural numbers. Then 
for all k £ J7* we define: 

■Kij(a;k) := / dxui(x; k)[(p a + k a )uj(x; k)] = {m(- ; k), (p a + k a )uj(- ;k)). (3-14) 
Jn 

Note that due to the phases presence in the eigenfunctions Uy's, we cannot be sure that the Ri,j's 
are continuous/diffcrcntiable at crossing points. But all these 'bad' phase factors will disappear 
when we take the traces (sec (|3.20|) and (|3.23j) below). 

We now can write the local traces of Proposition 13.21 in the following way: 

Proposition 3.3. Let (3 > and po > be fixed. Let fx^ = fJ-oo(P-,Po) € K be the unique solution 
of the equation Poo(P, e , 0) = pq. Then both quantities (|3.8|) and (|3.9[) can be rewritten as: 

1 1 f 

Tr L2(K 3){xaWcx,,i(AMoo,0)} = -7775^ 2J / dkC ^ >*> ^^O) 

1 1 j 1 ,...,j i =i Jn * 



an 



d 



±f4t s , (3.15) 

2«r 7 r (j^ ( k ) _ e ) (£ . 2 ( k) _ £) (£ . 3 ( k ) _ £) {e u (k) - £) 



Tr L2(R 3 ){x ^ 2(/3 ,^,0)}^4^{| i £dkJ- jfdf 



/r (%(k)-0 J 

where the functions Q* 3 k t— > ,j 2 -ja-j4. ( k ) an d ^* 3 k i— > C 31) j 2 (k) are defined by: 

(l;k)^ 4jl (2;k)} (3.17) 
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and 

CjuM :=7r Jlj2 (l;k)7r j2jl (l;k) + 7r J1 . j2 (2;k)7r j2 . J1 (2;k) = \n h J2 (1; k)| 2 + \% 3l . n (2; k) | 2 . (3.18) 

Wc do not give more details since this result is just a straightforward application of the following 
rather non-trivial technical lemma (recently proved in |16|): 

Lemma 3.4. Let j3 > and \i eK be fixed. For n, m <E N with m, ?i > 1, consider the local trace 
given by: 

j£L«n ■■= T ^(K 3 ){xo^da(/3,/, ; ^(^(o) -0-^(^(0) -er 1 ---Pa B (Hoo(o) 

TTiera under the assumption that V € C°°(T 3 ) we /lave: 
^i™La„ = E dk ^iJ2 (an; k) . . . 7t jnijl (a n ; k) 

f(p>;0 



ji,... j'„>i 



d£ LiLJ ^' s/ (3 19) 

It (E h (k) ~ m+1 (%(k) -€)•" (^n(k) - £) ' 
where all the above series are absolutely convergent and 7^(0!; k) is defined by (|3.14j) . 

3.3 Applying the residue calculus 

Consider the expression of the susceptibility at fixed density (|3.7j) in which the local traces are 
now given by (|3.15[) and (|3.16p . Remark that these quantities now are written in a convenient way 
in order to apply the residue theorem. Denote the integrands appearing in (|3.15[) and p,16[) by: 

1 (a r\ f(p\Mao)v) • r- M * 

(^(k)-0 2 (^ a (k)-0(^ 3 (k)-0(^(k)-0 
Note that Qj 1 ,j 2 {Pi f-oo'i • ) can have first order, third order, or even fourth order poles (in the case 
when ji = ]2). In the same way, ^)j 1< j 2 ,j3,0i(Pi ^°o'i 1 ) can have poles from the first order up to 
at most fifth order (in the case when j\ = ji = j'3 = j'4). Hence we expect that the integrals of 
fyh J2 j3 ,34 {Pi Moo ; • ) in (|3-15[) (rcsp. of Qj 1 J2 (/?, y.^ ; • ) in (|3.16p ) to make appear partial derivatives 
of f(/3,/i oc ; • ) with order at most 4 (resp. with order at most 3). But we will see below that the 
factor multiplying (<9|f)(/3, /Utx,; • ) is identically zero. 

Getting back to the susceptibility formula in p.7[) and by virtue of the previous remarks, we 
expect to obtain an expansion of the orbital susceptibility of the type (|3 . 1 1) . The next two results 
identify the functions tj u i{- ) coming from (|3.15j) and p,16|) : 

Lemma 3.5. The quantity defined by p,15[) can be rewritten as: 



1 1 f 9 1 

Tr LHR3) { X nW 00 , 1 (j3,^)}=--—Yl / dk E ^7 O 3 ' ^ 5 E ^ ( k )) a n ,1 ( k ) ( 3 - 20 ) 



where for all ji £ N* and k G Q*, the functions Cij 1 ,3( - ) and ttj lt 2(-) are given by: 
Oh.aW-gjjl^dsk)! E^, (k) _ %(k) +1^(^)1 E^ (k) _^ (k) 



J2#jl .72 /jl 



* n-kW V 2jft( 7 r Jlj2 (2;k)7r 32jl (l;k)) | 

T 31J1 (UK,i 1 (2,k) X, ^(k)-^) J (3 - 21) 
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and 



/i\ _ _J_f fyi ,31 ,h <h (k) + Cj 1 j 2 j 1 j 3 (k) + Cj 1 ,j 2 j t (k) 

JlM j ' 2!\ 2- 2- (S, 2 (k) -^,(k))(^ 3 (k)-iUk)) 



J2=l 



(k)- ( k )l ,ooo. 



(E j2 (k)-^(k))' 

J2#jl 

Note that it is possible to identify in (|3 . 20[) all the functions Oj 1; j(-) for ji > 1 and I <G {1,0} 
since such a result is based only on identities provided by the residue theorem. However, the 
number of terms is large and we will not need their explicit expressions in order to prove our 
theorem. 

Now we treat the next term. 
Lemma 3.6. The quantity defined by (|3.16[) can be rewritten as: 

Tr i 2 (R3) { Xn W 00 , 2 (/3, Atoo )} = -— J2 / dkEfl^'^^iM) 6 *.^) (3.23) 

I I ji=1 Jo* /=0 t 

where for all integers j\ > 1 and all k £ il* we have: 

b Jll3 (k) := i{|* jlA (l;k)| 2 + (2; k)| 2 }, (3.24) 

h n/\ . 1 l^i^(l;k)| 2 + 1^,^(2; k)| 2 1 

b - 2(k) ■= "2 £ ^ a (k)-^(k) + 2' (3 - 25) 

3'2/j'i 

, , n /„ v y> Kji,j 3 (i;k)| + Kii,j2(2;k)| 2 _, ni , 

b il>a (k) := -(2 - a) ^ — — 3^ , s G {0, 1}. 

j 3 =i (Sj 2 (k) -%(k)) 

j2 #ii 



Thus our Lemmas 13.51 and 13.61 provide an expansion of the type announced in (|3.ip , where the 
coefficients are given by: 

c rfl ,i(k):=aj 1 , l (k) + I» il ,i(k), Ze {0,1,2,3}. (3.26) 

In particular, for all integer j% > 1 and for all k 6 ST, the functions Cj lt 3(-) and are 
respectively given by: 



c ii,3 



M :-i{|* tA (l ik )|.(l + t J^egi 



J2 = 

h ^ji 



2 

(k) 



% M-kUr ^ V ^(-:,,^k;-.„ f: ;l:k;n ,„ 
"^(l-k^^k) ^(kj-i^k) J (3 - 27) 
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and 



2\£ ^(kj-i^k) ^ (^(k)-^^ 



oo oo 



EV^ Pjl ; jl , J2 , J3 (^) + C ji ,J 2 , jl , J3 (k) + £jl , j 2 , J3 , 31 (jjO I /o QO - ) 



J2=l J3 = l 
J2#j'l 33&1 



where for all integers ji, j 2 , j 3 , j 4 eN', !!' 3k4 Cji^Ws^O*) i s defined in (|3.17l) . 

In order to conclude the proof of Theorem 13.11 it remains to use this last result (its proof is in 
the appendix of this section) : 

Lemma 3.7. For all integers ji > 1 and I € {0,1,2,3}, the maps f2* 3 k t— > a 7 - ll j(k) and 
n* 9 k i— > b Jli i(k) are bounded and continuous on any compact subset ofQ* where Ej 1 is isolated 
from the rest of the spectrum. 

Thus for all integers j\ > 1 and k £ fT, the maps Cj lt i(- ) appearing in (|3.1[) might be singular 
on a set with zero Lebesgue measure where Ej 1 can touch the neighboring bands. However, the 
whole integrand in (|3.1[) is bounded and continuous on the whole Q* because it comes from two 
complex integrals ( p.l5p and (|3.16[1 ) which do not have local singularities in k. 

3.4 Appendix - Proofs of the intermediate results 

Here we prove Lemmas 13.51 13.61 and 13.71 
Proof of Lemma 13.51 Let SI* 9 k h-> Cj 1 J2 J3 J4 (k) be the complex- valued function appearing in 

Cji,j2,j s ,uQ*) ■= {ftji,j 3 (1; k)7Tj 2 .j 3 (2; k) — ftj l! j 2 (2;\i)Ttj 2! j 3 (l;\i)} 

x {%A(2;k)% A (l;k) -7r j3>j4 (l;k)7r i4>il (2;k)}. (3.29) 

Note that this function is identically zero for the following combinations of subscripts: 

h = .72 = 33 = ji, ji = 32 = 33 + ji, ji = 33 = ji ^ 32- (3.30) 

Therefore the expansion of f|3 . 1 5[) consists of partial derivatives of f (/?, fioo ; ■ ) of order at most equal 
to three. On the other hand, since the functions Cj 1 _j 1 _j 1 j i (- ) and Cji . j 2 .ji , ji (• ) are identically equal 
to zero (see (|3.30l0 . the quadruple summation in (|3.15[) is reduced to : 

1 \ f At fG8,A»oo;0 



jl,... J4 = l 



W Jt (E n (k) - tY(E h (k) - (E J3 (k) - £) {E U (k) - f) 

f(/3,Moo;0 



S £ C —" (k, (^) /r* - £ )*fe( k ) - f) 

j.3#ii 



.71=1 32=1 
J2#j 




r 3 (i^,(k)-£) 2 (£ j2 (k)-£) 3 



f(/3,Moo;0 



at most 2 equal 
subscripts 



r (k) - (^2 (k) - e) (-E* (k) - (k) - f ) 

(3.31) 
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By applying the residue theorem in the first term of the right hand side of (|3.31[) we get: 



X C Ji JiJWi ( k ) ( 2 - ) [ d £ n^tX^/J rt\ t\ ~ ( '-' •' :k: 
»=i V 7 Jr (^i ( k ) - (^A ( k ) - »=i 

33+31 33+31 

<9 3 f 3 1 <9 2 f 

(0 i H oo ;E h QL)) + -- — ' ( J 9, /loo; £. i (k)) + 



3!£,,(k)-i^k)^3V- — ' ^ " 3!(^. 3(k) _^ i(k) ) 2 ae 

0* 



others terms involving -^r(/3, Moo; •); with Z < 1 (3.32) 



The function Cj 1 j 1 .j 3 .j 1 (• ) appearing in front of tL| ^ooj -E^ (k)) in (|3.32l) corresponds to J1) 3(- ) 



as 

since: 



Vk G ft*, C J1Jlj3jl (k) = |^ 1J1 (l;k)^ lj3 (2;k) - ix h , h (2; k)** , h (1; k) | . 

Note that the function Cj 1 j 1 j 3 j 1 (-) contributes to the term a^Qi too. 

By applying once again the residue theorem in the second term of the right hand side of (|3.31l) 
we obtain: 



00 f 1 1 <9 2 f 

= £ ( '— Jki ( - 2! (%(k) _i;. 2( k)) 2 ^ (AMoo;£;j2(k)) 



J2 = l 



32 

h+3i 



others terms involving — j^,/^; •), with Z < 1 p> . (3.33) 



The function Cj 1 .j 2 .j 2 .j 2 (- ) appearing in front of Jj| (/3, jU^; £j 2 (k)) contributes to ft^^G )■ 

It remains to isolate in (|3.3ip (where at most two subscripts are equal) all combinations which 
provide a second order derivative of f(/3,/i oc ; • ). These combinations are: 

ii = .72 7^ i3 , j'4 ; ii = .73 7^ J2 , j'4 ; ii = h ^ h , J3 ■ 

Finally, we once again apply the residue theorem and gathering all terms proportional with 
g|i(/3, Hoc] • )• The proof is over. ■ 
Proof of Lemma \3. 6\ By separating the cases ji = j'2 and ji 7^ j'2 , the double summation in 



the right hand side of (|3.16[) reads as: 



jl = lj2 = l 



,2i»/Jr (£„(k)-|) J (£ A (k)-{) 



?, C -«(i)/ r d f 



1 A f ^ fC8,Moo;0 



ji=i 

00 00 



r '(^(k)-C) 4 



f(/3,Moo;0 



y y c.-^ck)! -M / d^ ,v "'^" y . (3.34) 



3 

3l+3 



_(fc)(k)-e) J (E J2 (k)-e) 

By using the residue theorem in the first term of the r.h.s. of (|3.34l) : 
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This is only the one term which provides a third-order partial derivative of f(/3, ■ ). The rest 
of the proof is just plain computation using the residue theorem. We do not give further details. 
■ 

Proof of Lemma \3. 7| Let p a := —id a be the a component of the momentum operator with 
periodic boundary conditions in L 2 (Vt), a G {1,2,3}. Now assume that Ej 1 {k) is isolated and 
non-degenerate if k belongs to some compact K C O*. We have to investigate integrals of the 
type 

Tr i2(n) ^d£f(/3, M oo;0(/>(k) - £)-WMk) - C)" 1 ■ ■ WOO - (3-35) 

Let ko € X, and let Ti be a simple, positively oriented path surrounding Ej 1 (ko) but no other 
eigenvalue of h(ko). If |k— ko| is small enough, Li will still only contain Ej 1 (k). The projection 
n(k) corresponding to Ej 1 (k) is given by a Ricsz integral. We have: 

H(k) = hjv teihM-z)- 1 , (3-36) 
and is continuous at ko in the trace norm topology. Moreover, 

n(k)(Mk) - er 1 = F ,\ M k), (i - n(k))(&(k) - zy 1 = ^-.f dz—(h(k) z)-\ 

E h (k) - £ 2m J Fl z-£ 

(3.37) 

Clearly, H(k)(h(k) — £) _1 is analytic in £ in the exterior of Li. We can decompose the integral on 
r in (|3.35[) as a sum of three integrals, one of which being on a simple contour T2 around Ej 1 (ko), 
completely surrounded by T±. The other two integrals will never have Ej 1 (k) as a singularity, so 
they cannot contribute to the formula of Oj lt z(k). On the other hand, in the integral on T2 we can 
replace the resolvents with the decomposition in Q3.37P and use the fact that (1 — II(k))(/i(k) — £) _1 
is analytic if £ lies inside Y2- Now one can apply the Cauchy residue formula. For example, we 
can compute the integral in which we have II(k) at the extremities, and (1 — n(k)) in the interior; 
in that case Ej 1 = Ej 1 (k) will be a double pole: 

Tr L2(fi) f d£f(£)n(k)(/i(k) - - - n(k)) ■ ■ -p^(Mk) - e) _1 n(k) 

= 2 7 u{(%f)(£ J1 (k))Tr i2(n) (n(k)p ai (ft(k) - EjJ-^l - n(k)) • • -p Q4 n(k)} + (3.38) 

+ f(E h ( k ))^Tr i2(n) {u(k) Pai (h(k) - o-^i - n(k)) • • • (ft(k) - o _1 (i - n(k)) PQ4 } ?=£ . i(k) }■ 

Thus one contribution to o^^k) will be: 

Tr i2(n) {n(k)p Q1 (h(k) - ^j-^i - n(k)) ■ • • Pct4 n(k)} . 

This expression does not use eigenvectors, only resolvents and projectors. Since Ej 1 is continuous 
at ko, the map 

k^(i-n(k))(Mk)- J B 31 (k))- 1 = -L [ dz — J—^)-*)- 1 

2tti J Fl z - (k) 

is operator norm continuous at k , and the map k i-> n(k) is continuous in the trace norm. By 
using standard perturbation theory (see e.g. [IS]), the same holds for the maps: 

k^ (l-n(k))(/i(k) -^(k))- 1 ^, and kH+p Ql (l-n(k))(/i(k) -^(k))-^. 

Thus the trace defines a continuous function; all other coefficients can be treated in a similar way. 
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4 The zero-field susceptibility at fixed density and zero 
temperature 

In this section, we separately investigate the semiconducting and metallic cases from the ex- 
pansion (|3.1I) . In particular, we prove Theorem ll.2l (i) and (ii). 

4.1 The semiconducting case (SC)- Proof of Theorem 11.21 (i) 

By using that ^fd{P, Mi = — /3~ 1 c'^f(/3, ^; (|3.1[) can be rewritten as: 



eVl 1 



cj 2 (2tt) 3 

E / ^{E% £ (^^~5^( k )) c i I .^( k )-|K^M<»;^ 1 (k))c A ,o(k)}. (4.1) 

3 - 1= l ,/n * ^ i=0 ^ P ' 

From (|4.ip . the proof of Theorem 11.21 (i) is based on two main ingredients. The first one is that 
for any fixed p > Eo we have the following pointwise convergences: 

lim -f(/3,fj,;Q = (ft - Qx[e ,h](€), Jim fFD(P, M5 = X[e ,A (0> v C G [#0, oo) \ {//}, 

p— J-oo p p— )-oo 

(4.2) 

while in the distributional sense: 



lim r^S,^) = liffl^ii;^^^). (4.3) 

The second ingredient is related to the decay of the derivatives of the Fermi-Dirac distribution: 
for all d > and for all jgN*, there exists a constant Cj^ > such that 



sup 

|£-/i|>d>0 



0f» 



• C^e ' • (4-4) 



Now assume that we are in the semiconducting case with a non-trivial gap, that is there 
exists N € N* such that lim^oo /Zoo(/3, po) = (max£jv + min£jv+i)/2 = £f{po) and max£/v < 
minfjv_|_i. Since the Fermi energy lies inside a gap, all terms containing derivatives of the Fermi- 
Dirac distribution will converge to zero in the limit j3 — > oo. Here (|4.4[) plays a double important 
role: first, it makes the series in j\ convergent, and second, it provides an exponential decay to 
zero. Then by taking into account (|4.2p . wc immediately get (|1.15[) from (|4.1j) . ■ 



4.2 The metallic case (M)- Proof of Theorem 11.21 (ii) 

Now we are interested in the metallic case. The limit j3 — > oo is not so simple as in the previous 
case, because the Fermi energy lies in the spectrum. The starting point is the same formula (|3.1j) . 
but we have to modify it by getting rid of the third order partial derivatives of f in order to make 
appear a Landau-Peierls type contribution. However, this operation needs the already announced 
additional assumption of non-degeneracy (which will provide regularity in k) in a neighborhood 
of the Fermi surface: 

Assumption 4.1. We assume that there exists a unique JVeN" such that lim^oo Pooifi, Po) = 
£f(po) S (minf jvj maxf jv) 7 which means that the Fermi energy lies inside the Nth Bloch band 
Em- We also assume that the Fermi surface defined by Sf '■= {k G 0* : i£jv(k) = £f (po)} is 
smooth and non- degenerate. 
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Recall that En (k) is supposed to be non degenerate outside a (possibly empty) zero Lebesgue 
measure set of k-points. Our assumption leads to the following consequence: 



dist{£ F (Po), Uf^ 1 ^} = di > 0, dist{£ F ( Po ), ^>T=N+i £ j} = d 2 > 0. (4.5) 

Note that the minimum of the lowest Bloch band E\ is always simple. If the density po is small 
enough then Assumption 14 . 1 1 is automatically satisfied since the Bloch energy function k h-> £a(k) 
is non-degenerate in a neighborhood of k = (see e.g. [28]). 

In fact, the non-degeneracy assumption is indispensable for to use of the regular perturbation 
theory in order to express the functions defined by (|3.24[) , (|3.25[) and (|3.21[) (only in the case where 
ji = N) with the help of the partial derivatives of E^(-) with respect to the fc^-variablcs, for k in 
a neighborhood of the Fermi surface: 



= n N , N (i;k), ie {1,2,3}, (4.6) 



dE N (k) 
dh 

d 2 E N (k) 1 , ~ |7r ijJV (i;k) 12 



>^ ,.X.^ ^{1,2,3}, (4.7) 



dkf ^EnM-ew 



d 2 E N {k) = ~ 2Sft{7r J -, JV (l;k)7r Af , J -(2;k)} = d 2 E N (k) 
dkxdk2 ^ £?jv(k)-£j(k) <9fc 2 <%i ' 1 ' j 



Such identities have been studied in [IB]. Note that the above series are absolutely convergent if 
the potential V is smooth enough (|16j). 

Now using Assumption l4.il we can group the coefficients corresponding to the third and second 
order derivatives of f appearing in (|3.1[) . This operation allows us to isolate a Landau-Peierls type 
contribution (the proof can be found in the appendix of this section): 

Proposition 4.2. Assume for simplicity that £n is a simple band. Let f2* 3 k i— > CAr. 2 (k) and 
fi* 9 k i-> cjv,3(k) i/ie functions respectively defined by p.28|) and p.27|) with j\ = N . Then: 



J m dk E^(^~; j M k )) Cjv - i ( k ) 



d 2 f, a _ , j i a 2 Sjv(k) a 2 ^jv(k) i /a 2 Sjv(k)\ 

dk^(/3,^;^(k))|-^^^^--^^U) + a„, 2 (k)^ (4.9) 
where f2* 9 k n- aj 1:2 (k) are i/ie functions defined in (|3.22D . 



From (|3.ip and Proposition 14.21 we get an expansion for the orbital susceptibility at fixed 
density po > and inverse of temperature (3 > 0: 

Proposition 4.3. Assume for simplicity that £n is a simple band. For every ji G N* there 
exist four families of functions Cj-i,i{-) with I £ {0, 1, 2, 3}, defined on Q* outside a set of Lebesgue 
measure zero, such that the second integrand below is bounded and continuous on Q* : 



X(p,po) 



dk 



1 



1 



12(3 (2tt) 3 



dk0(/3,/i oo ;£; w (k)) 



d 2 E N (k)d 2 E N (k) f d 2 E N (k) 
dkidk 2 



dk 2 



dk 2 



h=i 1=2 ? 



i, —1 l— n ^ J J 



ii = l i=0 



(4.10) 
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where by convention (d®f)(j3) fioo] ') '■— f(fiiV>oo\') and: 



oo oo 



t (v\ ■ o„ fir\ — CN > N >h,h(k) +Cjv,j 2 ,jy,j 3 (k) +C Nij2 ,j 3 . N (k) 

T N (k) .= -2a„, 2 (k) - ^ ^ (^ (k )-^( k ))(^( k )-^( k) ) 



, Cj 2 ,N,N,N(k) -Cjy,jy,j 2| jy(k) ^ 

f£ x ~ (E j2 (k)-E N (k)) 2 



Note that we can use identities provided by the regular perturbation theory in order to express 
the functions Cj lt i (as well as J~n) appearing in (|4.10[) in terms of derivatives of Ej and Uj w.r.t. 
the k-variablc. But this formulation will only hold true outside a set of k-points of Lcbesguc 
measure zero, while the formulation involving Ttij's is more general, physically relevant, providing 
us with bounded and continuous coefficients on fi* (see Lemma I3.7j) . Finally keep in mind that 
the main goal is the Landau-Pcicrls formula, and it will turn out that only the factor multiplying 
the second partial derivative of f will contribute to it. 

In order to complete the proof of Theorem ll.2l (ii), it remains to take the limit when (3 — > oo 
in (|4.10[) . Since the Fermi energy lies inside the band £n and it is isolated from all other bands, 
then using (|4.3j) and (|4.4j) we have: 

lim - dk ^ £ ^T(/3^oo(/3,p );^(k))c,- ; (k) = 



3=1 1=2 



and 



L dk w^^ ] n{ )) \^ — dki — {-Wdkv) -^ N{k \ 

f d<r(k) f d 2 E N (k)8 2 E N (k) f d 2 E N (k) \ 2 \ 
J SF \VE N (k)\\ dkf dk 2 2 {dhdkzj 

where Sf denotes the Fermi surface. Using these two identities together with (|4.2|) in (|4. 10[) . we 
obtain (fl~T6|) . 

4.3 Appendix - Proof of Proposition 14.21 

Using (|3.26[) we get: 

/ dkV-i(ft^;£«(k))c w (k)= dk— I(/3, Aloo ;i? A r(k))a A r i2 (k) 
Jo,* l=2 °t, Jn* °S 

+ J a , dk £ 0(AMoo;^(k))b Wi/ (k) + 0(/3, A/oo ;£; A r(k))a A r i3 (k) 
Using ()4.6|) and (|4.7|) . the functions bjv,i(-)) ' € {2,3}, can be rewritten as: 

. 1 ffdE N Q&)\* , ^a^(k) 

bjv,3(kj = — 



3! \ V / V dfc 

'iv, 2 (k) - - 2? ( - 2 {-^kf- - 1 )~2 [~^kT~ 

n(d 2 E N (k) d 2 E N (k)° 



2!2\ 9fc2 
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Since En(- ) e C 2 (M 3 /(27rZ 3 )), a simple integration by parts gives us: 

Vie {1,2}, / dk t 

J — 7T 



dh dp 



dh 



dh^,^;E N (V)^V (412) 



whence: 



1 ^ dk (/?, Moo ; (k)) b iv.s (k) = ~ j ^ dk (/3, Moo ; Sjv (k)) 



dkl 



and: 



/• 3 chic ii f 

J nt dk J2d^(P>^EN (k)) f»JV,i (k) = gj 5 /„. dk ^ ^ ; Sjv (k) ) 



dfc 2 

d 2 E N (k) t ^^(k)) 



<9fc 2 



On the other hand, using (|4.6[) , (|4.7[) and (|4.8|) . the function 0jv.3( - ) can be rewritten as : 



3fc 2 2 j 
(4.13) 



Note that by a simple integration by parts: 

dE N (k) d 3 f 
dkj dp 

d 2 f d 
dk J 7^(P^oo;E N {k)) — 



J — 7T 



( dE N (k) \ ( dE N (k) \ f _ 8 2 E N (k) \ 
\ Oh J\ 8k 2 )\ dhdk 2 J 

dE N (k) d 2 E N (k) 



-)2 



(/3,/ioo;£;v(k))- 
dE N (k) d 2 E N (k) 



dkn 



dk 2 



(4.14) 



(4.15) 



dkj 



dk 2 



O T 



J L ^ J 

d 2 E N {k) d 2 E N (k) dE N {k) d d 2 E N (k)\ 



dk 2 



dk 2 



dkj dkj dk 2 



7 



£d*,||(A fe ; £ „( k )){^ 



d 2 E N (k) d 2 E N (k) dE N {k) d d 2 E N {k) 



dk 2 



dkj dki dkj dkj 



By virtue of P~T3|) . using P~T5j) and dHT5J), we get: 



J dk^|(/3, Moo ; J B JV (k))ai V ,3(k) = ii^dk0(/?, Aioo ;S JV (k))|2 



a 2 ^(k) a 2 £jv(k) 



<% 2 



<9fc 2 



dE N (k) d d 2 E N (k) dE N (k) d d 2 E N {k) d 2 E N (k) d 2 E N (k)\ 



dki dk 2 dk\dk 2 ' dk 2 dki dk 2 dk\ dk\ dk 2 J 

d 3 f( a „ n.^ 9E N (k) d E N (k) d 2 E N (k) 

dk 2 dk\dk 2 



(4.16) 



Finally, by a last integration by parts: 

dE N (k) d 2 f 
dk 3 dp 



Vi^je {1,2}, r 

J —TV 



dkj 



(fru^ENik)) 



d d 2 E N (k) 



dh dkj dkj 



dkj± 

oh 



dE N (k) d 2 f 



dkj 



dkj dp 

d 2 E N {k) d 2 f 



(/^/Zo^-E^k)) 



d 2 E N (k) 
dkj dki 



dkidkj dp 



(/3,Mcx=;-B JV (k)) 



dE N {k) dE N (k) d 3 f 
dkj dh dp 



(/3,/Zoo;£W(k)) 



d 2 E N (k) 
dkj dki 
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Then (|4.16l) is reduced to 
a f 



dk 



[p, /i oo ;£ , Af(k))0Ar i3 (k) = 



1 1 

3! 2 



dk 



(^,//oo;Sjv(k)) 



^^(^^^(k) J d 2 E N (k) \ 2 d 2 E N (k) d 2 E N (k) \ 
<9fc 2 ^ dhdh J dk 2 dk 2 J 



(4.17) 



By adding (|4~T3l) to (|4~T7)) we get : 

/ dk[vH(/3,/ioo;^(k))b Wi/ (k) 

d 2 f 



<9 3 f 

^3 (/3, Moo;-Ejv(k))ajv,3(k) 



1 

3! 



k^(/3,^oo;-B7v(k)) 



9 2 ^iv(k ) ^^(k) 



dk* 



d 2 E N (k) Y 
dk\dk2 J 



and we are done. Note that the proof does not work if En can touch other bands because we loose 
regularity. In that case the integration by parts have to be done across a tubular neighborhood 
of the Fermi surface Sp, the price being the apparition of some extra terms. These terms will 
though disappear in the limit (3 — > oo because they will decay exponentially with j3. ■ 



5 The Landau-Peierls formula 

The aim of this section is to establish an asymptotic expansion of (|1.16p in the limit of small 
densities (p —> 0). Here we prove the expansion (|1.17[) . of which (|1.18[) is a particular case which 
has allready been suggested by T. Kjeldaas and W. Kohn in 1957 |26j . 

5.1 Proof of Theorem 11.21 (Hi) 

Let us recall that E = min^o* Ei(k) = E 1 (0), and £a(k) is non degenerate near the origin 
with a positive definite Hessian matrix (see e.g. |28j). The same reference insures the existence of 
the following quadratic expansion of £i(k) for k —> 0: 

E 1 {k) = E Q + ^k T 

As the Hessian matrix is symmetric, then up to a change of coordinates this quadratic expansion 
can be rewritten as : 

1 3 ^2 

£i(k) =E + ~Y,^+O(k 4 ) whenk^O (5.1) 

1 i=l m i 

where \-Kr \ , ^ . arc the eigenvalues of the inverse effective-mass tensor. 

Consider the assumption of weak density pa £ (0, 1). In this case the Fermi energy defined by 
(|1.12j) lies in the interval (Eq, maxkgn* E'i(k)). When po — ^ it follows that £f(po) converges to 
Eq. The k-subset of f2* where Eq < £i(k) < £p(po) is therefore only localized near the origin. 

From (|5.ip we get the following asymptotic expansion of £f{po) — Eq when po — > (the proof 
is given in the appendix of this section): 

Proposition 5.1. When po — > ; we have the following expansion: 

£ F (p )-E = sp$+O p 3 , s: =L-J-( — - . (5.2) 
v ' 2 \m\m1 2 rn% i ) 

In the particular case when m* = m* > for i G {1,2,3} and by setting kp := (6ir 2 po)% : 

£ F (p )-E Q = -^k 2 F + O(k F ). (5.3) 



2 E 1 
dkjdk. 



(0) 



k + C(k 4 ) whenk^O 



. 1 <i.i<3 
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Before proving Theorem 11.21 (Hi), we need one more technical result (its proof is also in the 
appendix of this section): 

Lemma 5.2. Assume that Ei(k) remains non- degenerate on the ball _B Co (0) :— {k G CI* : |k| < 
eo} with eo > small enough. Consider any continuous function F : B eo (0) — > C. Then when 
Pq — > we have the following asymptotic expansions: 

f i^wuM F (V=M +o(pl) with A := y/^*^*AV2TrF(0)^ (5.4) 
Js F |V^i(k)| 

and 

f dk X[Eo , £F(po)] (E 1 (k))F(k) = Bp + o(po) with B:= >lm>^F(0) S i, (5.5) 
Jn* 6 

where s is the coefficient defined in ([57 



Now we are ready to prove the Landau-Pcicrls formula in (|1.17[) . For this, consider the formula 
(I1.16[) . Remember that -E'i(-) is non-degenerate and analytic in a neighborhood of the origin. Let 
us concentrate ourselves on the first term appearing in (|1.16[) : 

eVl 1 f da(k) [ ^(k) ffi^k) / 9 2 gx(k) \ 2 \ 
c) 12 (2tt) 3 J Sp |V^i(k)| I dkl dkl \dk1dk2) [a - K>) 

since only this term will have a nonzero contribution to the leading term in (|1.17[) . The other 
term will go to zero like po; this can be shown using (|5.5p . (|5.ip . and the fact that the coefficients 
Ci : i and C1.0 are continuous near (see Lemma f3.T|) . 
Now consider the following function: 

By taking into account that J-\(- ) = — 2ai.2(- ) (see (|4. 1 1[) ) and by virtue of Lemma [3 .71 F(-) is 
continuous near the origin. According to (|5.4p , the only thing we need to do is to compute F (0). 
The determinant of the Hessian matrix gives after a short computation: 

^(k)^(k) ^ l( k)v_ 1 whcnk ^ 



(67r) 3 p§ + o(pq ) when p 0. 



<% 2 <9fcf V dk\dk2 J m\m\ 
Thus we can write: 

2 1 



<*m(po) = -( - 
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3^1(0) 



The only thing we have left to do, is proving that .Fi(O) = 0. The definition of T\ can be found 
in (|4.1ip . while the coefficients entering in its definition are defined in Q3.29p . 
Let us start by showing that for all integers j2: J3 > 2 we have: 

Ci,i,hdM =Ci,,w(<>) = C ia ,i,i,i(0) =C lilj2jl (0) =0. 

Indeed, in the expression of each of these functions it is possible to identify a factor of the type 
•^1,1(0?; 0), a € {1, 2} which are nothing but partial derivatives of Ei at the origin, thus they must 
be zero. It follows that: 

Tx (0) = £ £ Wh (°) - ^wft (o) - ^ (o)) ■ (5 - 8) 



2G 



Since: 

Chh,hh(0) = *i,h (I! 0)7r i2>1 (2; 0)*i Ja (2; 0)7r is>1 (l; 0) + 7t hh (2; 0)7^(1; 0)^ lj3 (l; 0)7r i3>1 (2; 0) 
- 7T 1Ja (2; 0)7r jai i(l; 0)7T lij3 (2; 0)7r ia ,i(l; 0) - tt 1j2 (1; 0)7r ia , 1 (2; 0)7r lj3 (1; 0)7r J - 3jl (2; 0), 
then (15.81) can be rewritten as : 



•Fi(O) =2 



^ ^ (2; 0)^,1(1; 0) 2 _ 7r u (2; 0)rc 3 , x(l; 0) \ 2 _ ^ 7^(2; 0)^, r(l; 0) x : 



^(oj-^iCo) 



E,-(0)-Bi(0) 7 E j (0)-E 1 (0) 



(5.9) 

But for k = we may choose all our eigenfunctions ui(- ; 0) to be real. It means that for all 
integers j > 2 and a € {1,2}, the matrix elements 7i"i j(a;0) are purely imaginary. As a result, 
the sums in (|5.9p are real numbers and cancel each other, thus ^1(0) =0. ■ 

5.2 Appendix - Proofs of intermediate results 



Here we prove Proposition 15.11 and Lemma 15.21 

Proof of Proposition [5711 In (|5.1|) use the change of variables fcj := , fci = . with i e {1,2,3}. 

v m * 

This gives: 

£i(k) := Eiiy/mfk) = E + ^{~k 2 + ~k 2 + fcg} + G(k 4 ). 
In spherical coordinates: 

#l(r, 6, <f>)=E + ir 2 + C(r 4 ) when r -> 0. (5.10) 

We would like to express r as a function of £1, and 0. Clearly, the equation E\(r{0, <fi), 9, 0) = 
E + A has a unique solution r(0, <f>, A) if A > is small enough. This solution obeys a fixed point 
equation of the type: 

r(0, 4>, A) = V2A[1 + O{r 2 (0, 0, A))] (5.11) 

which leads to the estimate: 

r(0, 0, A) = V2A[l + O(A)] when A ^ 0. (5.12) 



We can finally determine A (thus the Fermi energy) as a function po- By setting Q* 
it follows from (12.61): 



(27T) 

Using spherical coordinates: 



Po = y dkx[jj ,£; +A](^i(k)j. 



P o = V ^ 2 3 / d0 / dflsin^ / drr 2 + / d 



(2tt) 3 



lr r 2 



2A 



This is the equation we have to solve in order to find A as a function of po- Then by standard 
fixed point arguments we arrive at the estimate l|5.2p and we are done. I 
Proof of Lemma \5.2\ We only prove (|5.4j) . the other estimate being similar. As before, we prefer 
the new variables h = where i € {1,2,3}. Denote by E'i(k) = E'i(k), by F(k) = F(k) and 



with f2* := , Then we can formally write: 



l^ F(k, = / . dk{(&ta) - El(k » F(k) 

= yjm\m* 2 ml [ dkS(£ F {p Q ) - £i(k))F(k) 



y/mlmtmt [ dg(k) F(k). (5.13) 
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The quadratic expansion ()5.1[) implies |V^^i(k)| = |k| [l + C(k 2 )] when k->0. Then: 

/ I^fSm F(k) = V^* 2 m*F (0) / da(k) Ik)" 1 [1 + 0(1)]. (5.14) 

Using spherical coordinates, let us denote as before by r(9,(j),po) the unique root of the equation 
Ei (r(8, <fi, po), 6, cp\ = £p(po). Then f|5 . 14[) can be rewritten as: 



2jv 



F(k) = y/mlm*m* 3 F(0) / d</>/ d6 sin(0) r(0, 0, p ) [l + o(l)] 



dcr(k) 
/ Sf |VSi(k)| 

Now by setting A := £p(p ) — E and by using (|5 . 1 2[) when A — > 0: 

/ I^f F ^ = V^W^l 4V2^F(0) VA [l + o(l)]. 
Finally, we use (|5.2I) and the proof is over. 
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